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Abstract 



, For a quantum system, a density matrix p that is not pure can arise, via 



averaging, from a distribution /x of its wave function, a normalized vector belonging 
Q ! to its Hilbert space J^. While p itself does not determine a unique p, additional 

I facts, such as that the system has come to thermal equilibrium, might. It is thus 

not unreasonable to ask, which p, if any, corresponds to a given thermodynamic 
0^] ensemble? To answer this question we construct, for any given density matrix 

p, a natural measure on the unit sphere in J^, denoted GAP{p). We do this 
I using a suitable projection of the Gaussian measure on J^f with covariance p. 

qh| We establish some nice properties of GAP{p) and show that this measure arises 

naturally when considering macroscopic systems. In particular, we argue that it 
' is the most appropriate choice for systems in thermal equilibrium, described by 

;h ■ the canonical ensemble density matrix = (1/Z) exp(— GAP(p) may also 



be relevant to quantum chaos and to the stochastic evolution of open quantum 
systems, where distributions on are often used. 

Key words: canonical ensemble in quantum theory; probability measures on 
Hilbert space; Gaussian measures; density matrices. 



*Departments of Mathematics and Physics, Hill Center, Rutgers, The State University of New Jersey, 
110 Frelinghuysen Road, Piscataway, NJ 08854-8019, USA. E-mail: oldstein@math.rutgers.edu 

^Departments of Mathematics and Physics, Hill Center, Rutgers, The State University of New Jersey, 
110 FreHnghuysen Road, Piscataway, NJ 08854-8019, USA. E-mail: lebowitz@math.rutgers.edu 

^Mathematisches Institut, Eberhard-Karls-Universitat, Auf der Morgenstelle 10, 72076 Tiibingen, 
Germany. E-mail: tumulka@everest.mathematik.uni-tuebingen.de 

^Dipartimento di Fisica dell'Universita di Gcnova and INFN sezione di Genova, Via Dodecaneso 33, 
16146 Genova, Italy. E-mail: zanghi@ge.infn.it 



1 



1 Introduction 



In classical mechanics, ensembles, such as the microcanonical and canonical ensembles, 
are represented by probability distributions on the phase space. In quantum mechan- 
ics, ensembles are usually represented by density matrices. It is natural to regard these 
density matrices as arising from probability distributions on the (normalized) wave func- 
tions associated with the thermodynamical ensembles, so that members of the ensemble 
are represented by a random state vector. There are, however, as is well known, many 
probability distributions which give rise to the same density matrix, and thus to the 
same predictions for experimental outcomes j2Sl sec. IV.3].^ Moreover, as emphasized 
by Landau and Lifshitz [ini sec. 1.5], the energy levels for macroscopic systems are so 
closely spaced (exponentially small in the number of particles in the system) that "the 
concept of stationary states [energy eigenfunctions] becomes in a certain sense unrealis- 
tic" because of the difficulty of preparing a system with such a sharp energy and keeping 
it isolated. Landau and Lifshitz are therefore wary of, and warn against, regarding the 
density matrix for such a system as arising solely from our lack of knowledge about the 
wave function of the system. We shall argue, however, that despite these caveats such 
distributions can be both useful and physically meaningful. In particular we describe 
here a novel probability distribution, to be associated with any thermal ensemble such 
as the canonical ensemble. 

While probabihty distributions on wave functions are natural objects of study in 
many contexts, from quantum chaos [3^1123 to open quantum systems our main 
motivation for considering them is to exploit the analogy between classical and quantum 
statistical mechanics [2111 EH ESI CH CHI EE] • This analogy suggests that some relevant 
classical reasonings can be transferred to quantum mechanics by formally replacing the 
classical phase space by the unit sphere ^(^) of the quantum system's Hilbert space 
J^. In particular, with a natural measure fJ^idip) on S^{J^) one can utilize the notion 
of typicality, i.e., consider properties of a system common to "almost all" members of 
an ensemble. This is a notion frequently used in equilibrium statistical mechanics, as 
in, e.g., Boltzmann's recognition that typical phase points on the energy surface of a 
macroscopic system are such that the empirical distribution of velocities is approximately 
Maxwellian. Once one has such a measure for quantum systems, one could attempt an 
analysis of the second law of thermodynamics in quantum mechanics along the lines of 

^This empirical equivalence should not too hastily be regarded as implying physical equivalence. 
Consider, for example, the two Schrodinger's cat states ^± = (\E' alive ± '^dead)/ V^. The measure that 
gives equal weight to these two states corresponds to the same density matrix as the one giving equal 
weight to ^aiive and \E'dead- However the physical situation corresponding to the former measure, a 
mixture of two grotesque superpositions, seems dramatically different from the one corresponding to 
the latter, a routine mixture. It is thus not easy to regard these two measures as physically equivalent. 
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Boltzmann's analysis of the second law in classical mechanics, involving an argument 
to the effect that the behavior described in the second law (such as entropy increase) 
occurs for typical states of an isolated macroscopic system, i.e. for the overwhelming 
majority of points on with respect to fi{dip). 

Probability distributions on wave functions of a composite system, with Hilbert 
space Jif, have in fact been used to establish the typical properties of the reduced 
density matrix of a subsystem arising from the wave function of the composite. For 
example. Page ^H] considers the uniform distribution on S^{Jif) for a finite-dimensional 
Hilbert space J^, in terms of which he shows that the von Neumann entropy of the 
reduced density matrix is typically nearly maximal under appropriate conditions on the 
dimensions of the relevant Hilbert spaces. 

Given a probability distribution fi on the unit sphere S^{Jif) of the Hilbert space 
there is always an associated density matrix |2S]: it is the density matrix of the 
mixture, or the statistical ensemble of systems, defined by the distribution fi, given by 



For any projection operator P, tr {Pfj,P) is the probability of obtaining in an experiment 
a result corresponding to P for a system with a /^-distributed wave function. It is 
evident from ((H) that is the second moment, or covariance matrix, of p, provided p 
has mean (which may, and will, be assumed without loss of generality since ip and —ijj 
are equivalent physically). 

While a probability measure p on S^{J^) determines a unique density matrix p on 
via ((H), the converse is not true: the association p p^ given by ((Tj) is many-to- 
one. ^ There is furthermore no unique "physically correct" choice of p for a given p since 
for any p corresponding to p one could, in principle, prepare an ensemble of systems 
with wave functions distributed according to this p. However, while p itself need not 
determine a unique probability measure, additional facts about a system, such as that 
it has come to thermal equilibrium, might. It is thus not unreasonable to ask: which 
measure on J^{Jif) corresponds to a given thermodynamic ensemble? 

Let us start with the microcanonical ensemble, corresponding to the energy interval 
[E, E + 6], where 6 is small on the macroscopic scale but large enough for the interval 

^For example, in a fc-dimensional Hilbert space the uniform probability distribution u = u yi^^) over 
the unit sphere has density matrix pu — \l with / the identity operator on at the same time, for 
every orthonormal basis of the uniform distribution over the basis (which is a measure concentrated 
on just k points) has the same density matrix, p — jl. An exceptional case is the density matrix 
corresponding to a pure state, p = \ijj){ip\, as the measure p, with this density matrix is almost unique: 
it must be concentrated on the ray through -ijj, and thus the only non- uniqueness corresponds to the 
distribution of the phase. 
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to contain many eigenvalues. To this there is associated the spectral subspace M'e,&i 
the span of the eigenstates \n) of the Hamiltonian H corresponding to eigenvalues En 
between E and E + 5. Since J^e,s is finite dimensional, one can form the microcanonical 
density matrix 



with P^'j^ g = l[E,E+s]{H) the projection to J^e,s- This density matrix is diagonal in the 
energy representation and gives equal weight to all energy eigenstates in the interval 



But what is the corresponding microcanonical measure'! The most plausible answer, 
given long ago by Schrodinger |2(H l2T] and Bloch [2^], is the (normalized) uniform 
measure ue,s = g) on the unit sphere in this subspace. pe,5 is associated with ue,5 

via dH). 

Note that a wave function chosen at random from this distribution is almost cer- 
tainly a nontrivial superposition of the eigenstates \n) with random coefficients (n|\E') 
that are identically distributed, but not independent. The measure ue,5 is clearly sta- 
tionary, i.e., invariant under the unitary time evolution generated by H, and it is as 
spread out as it could be over the set ^{J^e,5) of allowed wave functions. This measure 
provides us with a notion of a "typical wave function" from J^e,5 which is very different 
from the one arising from the measure he,s that, when H is nondegenerate, gives equal 
probability {dim J^e,5)~^ to every eigenstate \n) with eigenvalue En G [E,E + 6]. The 
measure fiE,5, which is concentrated on these eigenstates, is, however, less robust to 
small perturbations in H than is the smoother measure ue,s- 

Our proposal for the canonical ensemble is in the spirit of the uniform microcanonical 
measure ue,s and reduces to it in the appropriate cases. It is based on a mathematically 
natural family of probability measures /i on ^{Jif). For every density matrix p on 
Jif, there is a unique member /i of this family, satisfying (0) for = p, namely the 
Gaussian adjusted projected measure GAP{p), constructed roughly as follows: Eq. 
(i.e., the fact that p^ is the covariance of p) suggests that we start by considering the 
Gaussian measure G{p) with covariance p (and mean 0), which could, in finitely many 
dimensions, be expressed by G{p){dip) oc eyiY>{—{ip\p~^\ip)) dip (where dip is the obvious 
Lebesgue measure on M').'^ This is not adequate, however, since the measure that we 
seek must live on the sphere ^(J^) whereas G{p) is spread out over all of J^. We 
thus adjust and then project G{p) to in the manner described in Section 121 

in order to obtain the measure GAP{p), having the prescribed covariance p as well as 

•^Berry |3] has conjectured, and for some cases proven, that such measures describe interesting 
universal properties of chaotic energy eigenfunctions in the semiclassical regime, see also ^JESl- It 
is perhaps worth considering the possibility that the GAP measures described here provide somewhat 
better candidates for this purpose. 




(2) 



[E,E + 5]. 
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other desirable properties. 

It is our contention that a quantum system in thermal equilibrium at inverse tem- 
perature (3 should he described by a random state vector whose distribution is given by 
the measure GAP{p[^) associated with the density matrix for the canonical ensemble, 

Pfi = P^,H,i3 = exp(-/3if) with Z := tr exp{~(3H). (3) 

In order to convey the significance of GAP{p) as well as the plausibility of our 
proposal that GAP{pp) describes thermal equilibrium, we recall that a system described 
by a canonical ensemble is usually regarded as a subsystem of a larger system. It is 
therefore important to consider the notion of the distribution of the wave function of a 
subsystem. Consider a composite system in a pure state ip G J^i ® and ask what 
might be meant by the wave function of the subsystem with Hilbert space J^i. For 
this we propose the following. Let {|g2)} be a (generalized) orthonormal basis of 
(playing the role, say, of the eigenbasis of the position representation). For each choice 
of 1^2), the (partial) scalar product {q2\'ip), taken in is a vector belonging to Mi. 
Regarding \q2) as random, we are led to consider the random vector \I'i G J^i given by 

^i=U{Q2\i^) (4) 

where N = N'{iIj,Q2) = \\{Q2\'ip)\\ ^ is the normalizing factor and \Q2) is a random 
element of the basis {lg'2)}, chosen with the quantum distribution 

P(g2 = g2) = ||(g2|^)ir. (5) 

We refer to as the conditional wave function [Hj of system 1. Note that becomes 
doubly random when we start with a random wave function in ® instead of a 
fixed one. 

The distribution of \E'i corresponding to (@|) and (0) is given by the probability 
measure on y{J^i), 

/ii(#i) = P(^i G d^i) = 5^11 (g2|^)|r 5(^1 -Ar(^, 92) (g2|^)) #1, (6) 

where 6{tp — (p) dip denotes the "delta" measure concentrated at 0. While the density 
matrix p^^ associated with pi always equals the reduced density matrix p^^"^ of system 
1, given by 

pr' = tr2|V^)(^| = $^(g2|V^)(V'|g2), (7) 

92 

the measure pi itself usually depends on the choice of the basis {|g2)}- It turns out, 
nevertheless, as we point out in Section IKTT] that pi{dipi) is a universal function of pf^ 
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in the special case that system 2 is large and ip is typical (with respect to the uni- 
form distribution on all wave functions with the same reduced density matrix), namely 
GAP^pY'^). Thus GAP{p) has a distinguished, universal status among all probability 
measures on with density matrix p. 

To further support our claim that GAPi^pp) is the right measure for p/j, we shall 
regard, as is usually done, the system described by pp as coupled to a (very large) 
heat bath. The interaction between the heat bath and the system is assumed to be 
(in some suitable sense) negligible. We will argue that if the wave function of the 
combined "system plus bath" has microcanonical distribution ue,5i then the distribution 
of the conditional wave function of the (small) system is approximately GAP{pp)\ see 
Section H 

Indeed, a stronger statement is true. As we argue in Section 15. 2| even for a typical 
fixed microcanonical wave function ip of the composite, i.e., one typical for ue^s, the 
conditional wave function of the system, defined in Q, is then approximately GAP{pp)- 
distributed, for a typical basis {|g2)}- This is related to the fact that for a typical 
microcanonical wave function ip of the composite the reduced density matrix for the 
system is approximately [H Note that the analogous statement in classical 

mechanics would be wrong: for a fixed phase point ^ of the composite, be it typical or 
atypical, the phase point of the system could never be random, but rather would merely 
be the part of ^ belonging to the system. 

The remainder of this paper is organized as follows. In Section |21 we define the 
measure GAP{p) and obtain several ways of writing it. In Section El we describe some 
natural mathematical properties of these measures, and suggest that these properties 
uniquely characterize the measures. In Section HI we argue that GAP{pp) represents the 
canonical ensemble. In Section El we outline the proof that GAP{p) is the distribution 
of the conditional wave function for most wave functions in J^i ® with reduced 
density matrix p if system 2 is large, and show that GAP{p/^) is the typical distribution 
of the conditional wave function arising from a fixed microcanonical wave function of 
a system in contact with a heat bath. In Section IHl we discuss other measures that 
have been or might be considered as the thermal equilibrium distribution of the wave 
function. Finally, in Section [71 we compute explicitly the distribution of the coefficients 
of a Gy4P(p^)-distributed state vector in the simplest possible example, the two-level 
system. 

2 Definition of GAP(p) 

In this section, we define, for any given density matrix p on a (separable) Hilbert space 
J^, the Gaussian adjusted projected measure GAP{p) on ^(^). This definition makes 
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use of two auxiliary measures, G{p) and GA{p), defined as follows. 

G{p) is the Gaussian measure on with covariance matrix p (and mean 0). More 
explicitly, let {|?t,)} be an orthonormal basis of eigenvectors of p and Pn the corresponding 
eigenvalues, 

P = ^Pn\n){n\. (8) 

Such a basis exists because p has finite trace. Let Z„ be a sequence of independent 
complex-valued random variables having a (rotationally symmetric) Gaussian distribu- 
tion in C with mean and variance 

nZn\^=Pn (9) 

(where E means expectation), i.e., lie Zn and ImZ„ are independent real Gaussian 
variables with mean zero and variance p„/2. We define G{p) to be the distribution of 
the random vector 

^^:=5^Z„|n). (10) 

71. 

Note that \E'*^ is not normalized, i.e., it does not lie in In order that lie in 

at all, we need that the sequence Zn be square- summable, l-^nP < oo. That this 
is almost surely the case follows from the fact that ^Yln l^nl"^ is finite. In fact, 

= ^E|Z„|2 = J^Pn = tip = 1. (11) 

71 n n 

More generally, we observe that for any measure p on with (mean and) covariance 
given by the trace class operator C, 

we have that, for a random vector \I' with distribution p, E||^|p = trC. 

It also follows that '^'^ almost surely lies in the positive spectral subspace of p, the 
closed subspace spanned by those |n) with p„ 7^ 0, or, equivalently, the orthogonal 
complement of the kernel of p; we shall call this subspace support(p). Note further that, 
since G{p) is the Gaussian measure with covariance p, it does not depend (in the case 
of degenerate p) on the choice of the basis {\n)} among the eigenbases of p, but only on 
P- 

Since we want a measure on J^(J^) while G{p) is not concentrated on S^{J^) but 
rather is spread out, it would be natural to project G{p) to S^{J^). However, since 
projecting to ^^{J^) changes the covariance of a measure, as we will point out in detail 
in Section 13.11 we introduce an adjustment factor that exactly compensates for the 
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change of covariance due to projection. We thus define the adjusted Gaussian measure 
GA{p) on by 

GA{p){d^) = ll^f G{pm). (12) 

Since E||\I/*^|p = 1 by pij). GA{p) is a probabihty measure. 

Let '^'^^ be a GA(p)-distributed random vector. We define GAP{p) to be the dis- 
tribution of 

with P the projection to the unit sphere (i.e., the normahzation of a vector), 

P:jr\{o}^y{j^), p(^) = ii^irV- (14) 

Putting (O differently, for a subset B C ^(^), 

G'AP(p)(P)=GA(p)(M+i?)= 1 G(p)(#)||^f (15) 



where M+P denotes the cone through P. More succinctly, 

GAP{p) = P, (G'A(p)) = GA{p) o p-i . (16) 

where P* denotes the action of P on measures. 

More generally, one can define for any measure p on Jif the "adjust-and-project" 
procedure: let A{p) be the adjusted measure A{p){dilj) = \\ip\\^ p{dip)] then the adjusted- 
and-projected measure is P*(y4(p)) = A{p) o P^^, thus defining a mapping P* o A from 
the measures on Jif with J p{dip) ||^/'||^ = 1 to the probability measures on We 
then have that GAP{p) = P,{A{G{p))). 

We remark that \1>'^^-^, too, lies in support(p) almost surely, and that P(\l/'^) does 
not have distribution GAP{p) — nor covariance p (see Sect. IH.lj ). 

We can be more explicit in the case that p has finite rank k = dim support (p), e.g. for 
finite-dimensional J^: then there exists a Lebesgue volume measure A on support (p) = 
C^, and we can specify the densities of G{p) and GA{p), 



dGA{p) 



2 



W = JZ . exp(-(^|p;»), (17b) 



dX tt'^ det p_| 

with p_|_ the restriction of p to support (p). Similarly, we can express GAP{p) relative to 
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the {2k — 1) -dimensional surface measure u on ^ (support (p)), 

oo 

= -1-^ [ rfrr^'=-ir^exp(-r^(^|p;V)) = (18a) 
au tt'^ det J 





P+'li^)-'-'- (18b) 



2n'' det p+ 
We note that 

GAP{pE,s)=UE,s, (19) 

where pes is the microcanonical density matrix given in Q and Ues is the microcanonical 
measure. 



3 Properties of GAP{p) 

In this section we prove the following properties of GAP{p): 

Property 1 The density matrix associated with GAP{p) in the sense of is p, i.e., 
Pgap{p) = P- 

Property 2 The association p GAP{p) is covariant; For any unitary operator U on 

U,GAP{p) = GAPiUpU*) (20) 

where U* = U^^ is the adjoint ofU and is the action ofU on measures, U^^ = fioU^^ . 
In particular, GAP{p) is stationary under any unitary evolution that preserves p. 

Property 3 If E Mx ® has distribution GAP{pi (g) P2) then, for any basis {1^2)} 
of ^2, the conditional wave function has distribution GAP{pi). ("GAP of a product 
density matrix has GAP marginal. ") 

We will refer to the property expressed in Property 3 by saying that the family of 
GAP measures is hereditary. We note that when G Jifi®J^2 has distribution GAP{p) 
and p is not a tensor product, the distribution of \l'i need not be GAP{p'^^'^) (as we will 
show after the proof of Property 3). 

Before establishing these properties let us formulate what they say about our can- 
didate GAP{pp) for the canonical distribution. As a consequence of Property 1, the 
density matrix arising from p, = GAP{pp) in the sense of (P) is the density matrix pp. 
As a consequence of Property 2, GAP{pp) is stationary, i.e., invariant under the unitary 
time evolution generated by H . As a consequence of Property 3, if \E' G ^ = M\ ® 
has distribution GAP^p^^u^p) and systems 1 and 2 are decoupled, H = Hi®l2 + Ii®H2, 
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where /j is the identity on then the conditional wave function of system 1 has 
a distribution (in J^) of the same kind with the same inverse temperature f3, namely 
GAP{pj^-^^Hi,/3)- This fits well with our claim that GAPi^pp) is the thermal equilibrium 
distribution since one would expect that if a system is in thermal equilibrium at inverse 
temperature (3 then so are its subsystems. 

We conjecture that the family of GAP measures is the only family of measures 
satisfying Properties 1-3. This conjecture is formulated in detail, and established for 
suitably continuous families of measures, in Section 

The following lemma, proven in Section ESI is convenient for showing that a random 
wave function is GAP-distributed: 

Lemma 1 Let Q be a measurable space, p a probability measure on Q, and \1' : i7 — > 
a Hilbert-space-valued function. //\&(co') is G{p)- distributed with respect to p{duj), then 
^{uj)/\\^!{uj)\\ is GAP{p) -distributed with respect to \\^! {uj)\\'^ p{duj) . 



3.1 The Density Matrix 

In this subsection we establish Property 1. We then add a remark on the covariance 
matrix. 

Proof of Property 1 . From we find that 

PGAPi,)= I GAPip){d,p)\^){^\=E(^\^^^^){^^^^ 

^ e( \\^^^\\-^ \^^^){^^'^\] = [ G'A(p)(#) ||V^||~2|^)(V;| 



^ J G(p)(#)|^)(^|=p 

because j^,G{p){dip) is the covariance matrix of G'(p), which is p. (A number 

above an equal sign refers to the equation used to obtain the equality.) □ 

Remark on the covariance matrix. The equation pcAPip) = P can be understood as 
expressing that GAP{p) and G{p) have the same covariance. For a probability measure 
/i on Jif with mean that need not be concentrated on ^{J^), the covariance matrix 
is given by 

C,= [ M#)I0)(0I- (21) 
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Suppose we want to obtain from fi a probability measure on y{J^) having the same 
covariance. The projection P^,n of fi to S^{J^), defined by P^^{B) = fi{M.'^B) for 
B C ^(^), is not what we want, as it has covariance 



However, P^,(^A{n)j does the job: it has the same covariance. As a consequence, a nat- 
urally distinguished measure on ^{Jif) with given covariance is the Gaussian adjusted 
projected measure, the GAP measure, with the given covariance. 



We establish Property 2 and then discuss in more general terms under which conditions 
a measure on S^{J^) is stationary. 

Proof of Property 2. Under a unitary transformation U, a Gaussian measure with covari- 
ance matrix C transforms into one with covariance matrix UCU*. Since || [/■?/' |p = WipW"^, 
GA{C) transforms into GA{UCU*); that is, [/\Efg^ and "^ucv equal in distribution, 
and since Ht/^^f?"^!! = ||\E'f?^||, we have that U'^q^^ and are equal in distribution. 

In other words, GAP{C) transforms into GAPiUCU*), which is what we claimed in 



3.2.1 Stationarity 

In this subsection we discuss a criterion for stationarity under the evolution generated 
hy H = 'Y^^ En 1^) (^1- Consider the following property of a sequence of complex random 
variables Z^. 

The phases when they exist, are independent of the moduli |Z„| ^^^^ 

and of each other, and are uniformly distributed on = {e*^ : G M}. 

(The phase exists when Z„ ^ 0.) Condition (j^^ implies that the distribution 

of the random vector \E' = Ylm^nW) is stationary, since Znit) = exp(— ii?„t//?,)Z„(0). 
Note also that (|^^ implies that the distribution has mean 0. 

We show that the Z„ = {n\'if'^^^) have property (j22I)- To begin with, the Z„ = 
(nl^*^) obviously have this property since they are independent Gaussian variables. 
Since the density of GA{p) relative to G{p) is a function of the moduli alone, also the 
Zn = satisfy Finally, since the | (n|\E''^^'^) | are functions of the |(ra|\I''^^)| 

while the phases of the {n\'if'^^^) equal the phases of the (n|\l/'^'^), also the Zn = 
(^n\^GAP^ satisfy (f^. 




3.2 GAP{p) is Covariant 



(EOl). 



□ 
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We would like to add that is not merely a sufficient, but also almost a necessary 
condition (and morally a necessary condition) for stationarity. Since for any the 
moduli \Zn\ = |(n|\l')| are constants of the motion, the evolution of \E' takes place in the 
(possibly infinite-dimensional) torus 

{j2\Zn\e''"\n):0<9^<2n]^ J] S\ (23) 

contained in S^{J^). Independent uniform phases correspond to the uniform measure A 
on Yin is the only stationary measure if the motion on Yin is uniquely ergodic, 

and this is the case whenever the spectrum {En} of H is linearly independent over 
the rationals Q, i.e., when every finite linear combination 'Ylm'^nEn of eigenvalues with 
rational coefficients r„, not all of which vanish, is nonzero, see [21I21]- 

This is true of generic Hamiltonians, so that A is generically the unique stationary 
distribution on the torus. But even when the spectrum of H is linearly dependent, 
e.g. when there are degenerate eigenvalues, and thus further stationary measures on 
the torus exist, these further measures should not be relevant to thermal equilibrium 
measures, because of their instability against perturbations of H [TTl [T] . 

The stationary measure A on Hn corresponds, for given moduli or, equiva- 
lently, by setting \Zn\ = p(£'„)^/^ for a given probability measure p on the spectrum of 
if, to a stationary measure Ap on ^(^) that is concentrated on the embedded torus 
(I23|) . The measures Ap are (for generic H) the extremal stationary measures, i.e., the 
extremal elements of the convex set of stationary measures, of which all other stationary 
measures are mixtures. 

3.3 GAP Measures and Gaussian Measures 

Lemma n is more or less immediate from the definition of GAP{p). A more detailed 
proof looks like this: 

Proof of Lemma By assumption the distribution fi o \['~^ of with respect to fi is 
G{p). Thus for the distribution of \1/ with respect to p'{duj) = W"^ {uj)\\'^ fi^du) , we have 
n'o^^J-\d^) = \\^ffio^-\d^) = ll^f G(p)(#) = GA(p)(#). Thus, P(^(cj)) has 
distribution P^GA{p) = GAP{p). D 

3.4 Generalized Bases 

We have already remarked in the introduction that the orthonormal basis {lq'2)} of J^2, 
used in the definition of the conditional wave function, could be a generalized basis, such 
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"continuous" basis, for which it is 

l2 = 

instead of the "discrete" notation 

h = Y^ 1^2) (g2| 

we used in ©-(HI). 

We wish to elucidate this further. A generahzed orthonormal basis {\q2) '■ <i2 ^ Q2} 
indexed by the set Q2 is mathematically defined by a unitary isomorphism 
L'^{Q2, dq2), where dq2 denotes a measure on Q2. We can think of Q2 as the configuration 
space of system 2; as a typical example, system 2 may consist of particles in a box 
A C M^, so that its configuration space is Q2 = A^^ with dq2 the Lebesgue measure 
(which can be regarded as obtained by combining N2 copies of the volume measure on 
M^).^ The formal ket \q2) then means the delta function centered at 52; it is to be treated 
as (though strictly speaking it is not) an element of J^2- 

The definition of the conditional wave function \E'i then reads as follows: The vector 
■ip G Jifi ® can be regarded, using the isomorphism ^ L'^{Q2, dq2), as a function 
'ip '■ Q2 ^ '^1- Eq. (jH) is to be understood as meaning 

^^=Mi){Q2) (24) 

where 

M = M[^,Q2) = \m2)\~^ 

is the normalizing factor and Q2 is a random point in Q,2i chosen with the quantum 
distribution 

P(Q2Gc/g2) = ||^(g2)|rc^g2, (25) 

which is how (0) is to be understood in this setting. As ip is defined only up to changes 
on a null set in Q2, ^1 may not be defined for a particular (^2- Its distribution in 
J^, however, is defined unambiguously by ()24|1 . In the most familiar setting with 
= L\Qi,dqi), we have that (^(Q2))(gi) = ^{qi,Q2). 

In the following, we will allow generalized bases and use continuous instead of discrete 

notation, and set {Q2\i^) = 4'{Q2)- 

^In fact, in the original definition of the conditional wave function in [^1, 92 was supposed to be the 
configuration, corresponding to the positions of the particles belonging to system 2. For our purposes 
here, however, the physical meaning of the (72 is irrelevant, so that any generalized orthonormal basis 
of can be used. 



appropriate to write 
J dq2\q2){q2\ 
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3.5 Distribution of the Wave Function of a Subsystem 



Proof of Property 3. The proof is divided into four steps. 

Step 1. We can assume that ^ — P{^'^^) where is a G74(p)-distributed random 
vector in = JTi J^. We then have that *i = Pi((Q2|*)) = Pi(((52|*^^)) where 
Pi is the normahzation in M[, and where the distribution of Q2, given ^f*^^, is 

'^^^ and Q2 have a joint distribution given by the following measure v on M' x Q2: 

u{d^P X dq^) = \\{q2m?G{p){d^P)dq2. (26) 

Thus, what needs to be shown is that with respect to Pi((g2|'0)) is GAP{pi)- 
distributed. 

Step 2. If ^ E ® ,y^2 G{pi ® P2)- distributed and q2 G Q2 is fi.xed, then the 
random vector f{q2) (92!^) G ^ with f{q2) = (92|P2|q'2)^^^^ is G{pi) -distributed. This 
follows, more or less, from the fact that a subset of a set of jointly Gaussian random 
variables is also jointly Gaussian, together with the observation that the covariance of 
{q2\^) is 

/ G{pi {q2\i>) {'ilj\q2) = {q2\pi ® P2\q2) = Pi {q2\p2\q2) ■ 

More explicitly, pick an orthonormal basis {Ir^i)} of Ml consisting of eigenvectors of 
Pi with eigenvalues pn}, and note that the vectors \ni,n2) '■— \ni) (8) 1^2) form an or- 
thonormal basis of M — Mx ® M2 consisting of eigenvectors of p\ ® p2 with eigenvalues 
Pni,n2 —Pn}pnl- Sincc the random variables ^ni,n2 '■— {'>T'i,n2\^) are independent Gaus- 
sian random variables with mean zero and variances E|Z„^,n2p = Pni,n2-i so are their 
linear combinations 



Z{i)m ■■= (^1 1/(^2) ^fe)) = f{q2) ^{q2\n2) Z, 



m,n2 



"2 



with variances (because variances add when adding independent Gaussian random vari- 
ables) 

TUiv |2 t2i \Sr^\i I \|2Tn,|v |2 (l)Sn2 K^2|?T'2)Ppi2 (i) 
^Z{l)n^ = f [q2)lJ{q2\n2)\^Zn,,n2\ = / i i \ = Pnl ■ 

{q2\P2\q2) 

Thus 7(^2) (?2|^) is G(pi)-distributed, which completes step 2. 
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Step 3. If^E M{ ® ^2 is G{pi (8> p-z)- distributed and Q2 G Q2 is random with any 
distribution, then the random vector f{Q2) (Q2|^) is G{pi)- distributed. This is a trivial 
consequence of step 2. 

Step 4. Apply LemmaEas follows. Let Vt = J^xQ2, ^{to) = ^(^^,^2) = /(g2) {q2\^), 
and p{dip x dq2) = G{p){dip) {q2\p2\<i2) dq2 (which means that q2 and ip are independent). 
By step 3, the hypothesis of Lemma [T] (for p = pi) is satisfied, and thus Pi{^) = 
Pi{{q2\il>)) is G'y4P(pi)-distributed with respect to 

\\^{u;)fp{du;) = f{q2mq2m\'G{p){d^) {q2\p2\q2) dq2 = v{du) , 

where we have used that P{q2) = {q2\p2\q2)^^- But this is, according to step 1, what 
we needed to show. □ 

To verify the statement after Property 3, consider the density matrix p = |$)($| for 
a pure state $ of the form $ = y/p^'ipn ® (pn, where {ipn} and {0„} are respectively 
orthonormal bases for J^i and and the Pn are nonnegative with 'YlmPn = 1- Then 
a (jAP(p)-distributed random vector ^ coincides with $ up to a random phase, and 
so p^i'^ = J2nPn \'^n){i^n\- Choosiug for {192)} the basis {(pn}, the distribution of \E'i is 
not GAP{pY"^) but rather is concentrated on the eigenvectors of p^*^'^. When the p„ are 
pairwise-distinct this measure is the measure EIGipf^) we define in Section KTH 

4 Microcanonical Distribution for a Large System 
Implies the Distribution GAP^pj]) for a Subsystem 

In this section we use Property 3, i.e., the fact that GAP measures are hereditary, to show 
that GAP{pp) is the distribution of the conditional wave function of a system coupled 
to a heat bath when the wave function of the composite is distributed microcanonically, 
i.e., according to ue,5- 

Consider a system with Hilbert space J^i coupled to a heat bath with Hilbert space 
J^. Suppose the composite system has a random wave function G Jif = Jifi ® 
whose distribution is microcanonical, ue,5- Assume further that the coupling is negligibly 
small, so that we can write for the Hamiltonian 

H = Hi®l2 + h®H2, (27) 

and that the heat bath is large (so that the energy levels of H2 are very close). 

It is a well known fact that for macroscopic systems different equilibrium ensembles, 
for example the microcanonical and the canonical, give approximately the same answer 
for appropriate quantities. By this equivalence of ensembles [EI, we should have that 
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pE,5 ~ Pp for suitable (3 = P{E). Then, since GAP{p) depends continuously on p, we 
have that ue^s = GAP{pE^s) ~ GAP{pp). Thus we should have that the distribution 
of the conditional wave function \l'i of the system is approximately the same as would 
be obtained when \E' is G'74P(p^)-distributed. But since, by ^Tj\ . the canonical density 
matrix is then of the form 

PP = Pje,H,l3 = P.%i,Hi,l3 ® P.je2,H2,l3 , (28) 

we have by Property 3 that is approximately GAP{pjf^^ Hi,i3)-distrihuted, which is 
what we wanted to show. 

5 Typicality of GAP Measures 

The previous section concerns the distribution of the conditional wave function aris- 
ing from the microcanonical distribution of the wave function of the composite. It 
concerns, in other words, a random wave function of the composite. The result there is 
the analogue, on the level of measures on Hilbert space, of the well known result that if 
a microcanonical density matrix (j21) is assumed for the composite, the reduced density 
matrix p^^'^ of the system, defined as the partial trace tT2PE,5, is canonical if the heat 
bath is large [THj- 

As indicated in the introduction, a stronger statement about the canonical density 
matrix is in fact true, namely that for a fixed (nonrandom) wave function ip of the 
composite which is typical with respect to ue,5, Pi"^ ~ Pjifi,Hi,i3 when the heat bath is 
large (see [3121]; for a rigorous study of special cases of a similar question, see |22I)-^ 
This stronger statement will be used in Section 15.21 to show that a similar statement 
holds for the distribution of \E'i as well, namely that it is approximately GAP{p^^^H^^/^)- 
distributed for a typical fixed ip G J^e,s and basis {|g2)} of J^. But we must first 
consider the distribution of \E'i for a typical ip G J^. 

5.1 Typicality of GAP Measures for a Subsystem of a Large 
System 

In this section we argue that for a typical wave function of a big system the conditional 
wave function of a small subsystem is approximately GAP-distributed, first giving a 
precise formulation of this result and then sketching its proof. We give the detailed 
proof in [S]. 

^It is a consequence of the results in that when dimJ^ oo, the reduced density matrix be- 
comes proportional to the identity on J^i for typical wave functions relative to the uniform distribution 
on S^{J^) (corresponding to ue,5 ior E — and H — 0). 
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5.1.1 Statement of the Result 

Let M' = Jifi ® where Jifi and J^2 have respective dimensions k and m, with 
k < m < oo. For any given density matrix pi on Jifi, consider the set 

M{p,) = {^e : pr'(^) = Pi} , (29) 

where p^i'^{ip) = ^211^) {ipl is the reduced density matrix for the wave function ip. There 
is a natural notion of (normalized) uniform measure Up-^ on M{pi); we give its precise 
definition in Section IB. 1. HI 

We claim that for fixed k and large m, the distribution pf of the conditional wave 
function of system 1, defined by (jlj) and (0) for a basis {1^2)} of is close to 
GAP{pi) for the overwhelming majority, relative to Wp^, of vectors ip G with reduced 
density matrix pi. More precisely: 

For every e > and every bounded continuous function f : ^(^) M, 

Up,[ilje^{pi):\pt{f)-GAP{pi){f)\<e'j^l asm ^00, (30) 
regardless of how the basis {\q2)} is chosen. 
Here we use the notation 

p{f):= j pidi;)f{^). (31) 

5.1.2 Measure on J^f Versus Density Matrix 

It is important to resist the temptation to translate Up-^ into a density matrix in Jif. 
As mentioned in the introduction, to every probability measure p on ^{J^) there 
corresponds a density matrix p^ in Jif, given by (0), which contains all the empirically 
accessible information about an ensemble with distribution p. It may therefore seem 
a natural step to consider, instead of the measure p = Up-^, directly its density matrix 
p^ = ^pi®/2, where I2 is the identity on J^. But since our result concerns properties of 
most wave functions relative to p, it cannot be formulated in terms of the density matrix 
p^. In particular, the corresponding statement relative to another measure p' 7^ p on 
^{Jif) with the same density matrix p^/ = p^ could be false. Noting that p^ has a basis 
of eigenstates that are product vectors, we could, for example, take p' to be a measure 
concentrated on these eigenstates. For any such state ip, Pi is a delta-measure. 
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5.1.3 Outline of Proof 



The result follows, by ©, Lemmas and the continuity of P*oA, from the corresponding 
statement about the Gaussian measure G{pi) on Jifi with covariance pi: 

For every e > and every bounded continuous f : M\ M, 

n,,|^G^(pi):|/if(/)-G(pi)(/)|<£}^l asm^oo, (32) 

where fif is the distribution of y/m {Q2\4') G J^i (not normalized) with respect to the 
uniform distribution of Q2 G {1, . . . , m}. 

We sketch the proof of and give the definition of Up-^ . According to the Schmidt 
decomposition, every ip G can be written in the form 

V' = ^ Ci O 0i , (33) 

i 

where {xi} is an orthonormal basis of J^i, {(pi} an orthonormal system in J^2, and the 
Cj are coefficients which can be assumed real and nonnegative. From (jH^ one reads off 
the reduced density matrix of system 1, 

Pr = Y.^■\X^){X^\. (34) 
i 

As the reduced density matrix is given, p^"^ = pi, the orthonormal basis {xi} and the 
coefficients Cj are determined (when pi is nondegenerate) as the eigenvectors and the 
square-roots of the eigenvalues of pi, and, ^(pi) is in a natural one-to-one correspon- 
dence with the set 0NS{Jif2, k) of all orthonormal systems {(pi} in Jif2 of cardinality 
k. (If some of the eigenvalues of pi vanish, the one-to-one correspondence is with 
ONS{J^, k') where k' = dim support (pi).) The Haar measure on the unitary group of 
defines the uniform distribution on the set of orthonormal bases of J^2, of which 
the uniform distribution on ONS{J^, k) is a marginal, and thus defines the uniform 
distribution Up^ on .^(pi). (When pi is degenerate, Up^ does not depend upon how the 
eigenvectors Xi of pi are chosen.) 

The key idea for establishing from the Schmidt decomposition (jH^ is this: pf 
is the average of m delta measures with equal weights, pf = m^^ J2q2 ^ipi{q2)y located at 
the points 

k 

ipiiq2) =^Ci^/m{q2\(pi)x^■ (35) 
1=1 

Now regard ip as random with distribution Wp^; then the ipi{Q2) are m random vectors, 
and /if is their empirical distribution. If the mk coefficients (?2|0i) were independent 
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Gaussian (complex) random variables with (mean zero and) variance m~^, then the 
ipi{(l2) would be m independent drawings of a G(pi)-distributed random vector; by the 
weak law of large numbers, their empirical distribution would usually be close to G'(pi); 
in fact, the probability that (/) — G{pi){f)\ < e would converge to 1, as m — > oo. 

However, when {0j} is a random orthonormal system with uniform distribution as 
described above, the expansion coefficients {q2\4>i) in the decomposition of the (pi's 

<Pi = J2(^2\<l)i)\q2) (36) 

Q2 

will not be independent — since the (pi's must be orthogonal and since \\(pi\\ = 1. Nonethe- 
less, replacing the coefficients {q2\(pi) in flHUj) by independent Gaussian coefficients aj(g2) 
as described above, we obtain a system of vectors 

0: = E«^fe)ig2) (37) 

92 

that, in the limit m — > oo, form a uniformly distributed orthonormal system: ||0^|| 1 
(by the law of large numbers) and {(p'i\(p'j) —>■ for i j (since a pair of randomly chosen 
vectors in a high- dimensional Hilbert space will typically be almost orthogonal). This 
completes the proof. 

5.1.4 Reformulation 

While this result suggests that GAP{pp) is the distribution of the conditional wave 
function of a system coupled to a heat bath when the wave function of the composite is a 
typical fixed microcanonical wave function, belonging to ^e,s, it does not quite imply it. 
The reason for this is that J^e,& has measure with respect to the uniform distribution on 
J^, even when the latter is finite-dimensional. Nonetheless, there is a simple corollary, 
or reformulation, of the result that will allow us to cope with microcanonical wave 
functions. 

We have indicated that for our result the choice of basis {lq'2)} of does not 
matter. In fact, while /if, the distribution of the conditional wave function \E'i of system 
1, depends upon both ip G and the choice of basis {1^2)} of the distribution 
of /if itself, when ijj is -distributed, does not depend upon the choice of basis. This 
follows from the fact that for any unitary U on 

{U-\2\ip) = {q2\h<^Ui,) (38) 

(and the invariance of the Haar measure of the unitary group of under left multi- 
plication). It similarly follows from ()H8|1 that for fixed ip G J^, the distribution of /if 
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arising from the uniform distribution v of the basis {|5'2)}, in the set 0NB{^2) of all 
orthonormal bases of J^, is the same as the distribution of /if arising from the uniform 
distribution of with a fixed basis (and the fact that the Haar measure is invariant 
under U ^ U^^). We thus have the following corollary: 

Let ip G ^ and let pi = tT2\il'){'ip\ be the corresponding reduced density matrix for 
system 1. Then for a typical basis {\q2)} of J^, the conditional wave function \E'i of 
system 1 is approximately GAP{pi)- distributed when m is large: For every e > and 
every bounded continuous function f : ^{Jifi) —>■ M., 

u{^{\q2)}eONB{J^2):\ptif)-GAP{p^){f)\<E^ -^1 as dim (^2) oo . (39) 



5.2 Typicality of GAP{pp) for a Subsystem of a Large System 
in the Microcanonical Ensemble 

It is an immediate consequence of the result of Section 15.1.41 that for any fixed micro- 
canonical wave function ifj for a system coupled to a (large) heat bath, the conditional 
wave function of the system will be approximately GAP-distributed. When this is 
combined with the "canonical typicahty" described near the beginning of Sectional we 
obtain the following result: 

Consider a system with finite- dimensional Hilbert space M\ coupled to a heat hath with 
finite- dimensional Hilbert space J^. Suppose that the coupling is weak, so that we can 
write H = Hi®l2 + Ii®H2 on = J^i^J^2, o,nd that the heat hath is large, so that the 
eigenvalues of H2 are close. Then for any wave function ip that is typical relative to the 
microcanonical measure Ue^s, the distribution pf of the conditional wave function "^i, 
defined by (jlj) and (0) for a typical basis {1^2)} of the heat bath, is close to GAP{pp) 
for suitable (3 = 13 {E), where pp = p.%i,Hi,i3- In other words, in the thermodynamic 
limit, in which the volume V of the heat bath and dim{J^2) 90 to infinity and E/V = e 
is constant, we have that for all e,6 > 0, and for all bounded continuous functions 
f : R, 

UE,s X u[{i;,{\q2)}) E ^(Jf) x 0NB{.jr2) : \ptif) - GAP(p^)(/)| < e] ^ 1 (40) 
where P = /3(e). 

We note that if {|5'2)} were an energy eigenbasis rather than a typical basis, the 
result would be false. 
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6 Remarks 



6.1 Other Candidates for the Canonical Distribution 

We review in this section other distributions that have been, or may be, considered as 
possible candidates for the distribution of the wave function of a system from a canonical 
ensemble. 

6.1.1 A Distribution on the Eigenvectors 

One possibility, which goes back to von Neumann JSFj P- 329], is to consider ^{dip) as 
concentrated on the eigenvectors of p; we denote this distribution EIG{p) after the first 
letters of "eigenvector"; it is defined as follows. Suppose first that p is nondegenerate. 
To select an ii^JG' (p)-distributed vector, pick a unit eigenvector |n), so that p\n) = Pn\n), 
with probability p„ and randomize its phase. This definition can be extended in a natural 
way to degenerate p: 

EIG{p)= J2 Pdim^p M^(^i), (41) 
pespcc(p) 

where J^p denotes the eigenspace of p associated with eigenvalue p. The measure EIG{p) 
is concentrated on the set IJ^ J^p of eigenvectors of p, which for the canonical p = 
P^,H,p coincides with the set of eigenvectors of if; it is a mixture of the microcanonical 
distributions u yi^,^^) on the eigenspaces of H in the same way as in classical mechanics the 
canonical distribution on phase space is a mixture of the microcanonical distributions. 

Note that EIG{pe,s) = ue,5, and that in particular EIG{pe,s) is not, when H is 
nondegenerate, the uniform distribution fiE,5 on the energy eigenstates with energies in 
[E, E + 6], against which we have argued in the introduction. 

The distribution EIG{p) has the same properties as those of GAP{p) described in 
Properties 1-3, except when p is degenerate: 

The measures EIG{p) are such that (a) they have the right density matrix: Peig(p) = 
p; (b) they are covariant: U^EIG{p) = EIG{U pU*); (c) they are hereditary at nonde- 
generate p: when = Jifi ® and p is nondegenerate and uncorrected, p = pi® p2, 
then EIG{p) has marginal (i.e., distribution of the conditional wave function) EIG{pi). 

Proof, (a) and (b) are obvious. For (c) let, for i = 1,2, \ni) be a basis consisting of 
eigenvectors of pi with eigenvalues pn]- Note that the tensor products ® \n2) are 
eigenvectors of p with eigenvalues piiiPn-z j 

and by nondegeneracy all eigenvectors of p 
are of this form up to a phase factor. Since an ii^/G(p)-distributed random vector \& is 
almost surely an eigenvector of p, we have \E' = e*®|A'^i)|A^2) with random Ni, N2, and G. 
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The conditional wave function is, up to the phase, the eigenvector \Ni) of pi occurring 
as the first factor in The probabihty of obtaining A^^i = ni is PniPn2 = Pn}- 

In contrast, for a degenerate p = pi ® P2 the conditional wave function need not 
be i?/G'(pi)-distributed, as the following example shows. Suppose pi and p2 are non- 
degenerate but Pnlpn2 = PrniPrn2 for souie Ui ^ uii] then an i?/G'(p)-distributed \E', 
whenever it happens to be an eigenvector associated with eigenvalue Pnlpnl, is of the 
form c\ni)\n2) + c'|mi)|m2), almost surely with nonvanishing coefficients c and c 5 as a 
consequence, the conditional wave function is a multiple of c|ni)(Q2|^2)+c'|mi)((52|"^2), 
which is, for typical Q2 and unless \n2) and |m2) have disjoint supports, a nontrivial su- 
perposition of eigenvectors \ni), \mi) with different eigenvalues — and thus cannot arise 
from the EIG{pi) distribution.^ 

Note also that EIG{p) is discontinuous as a function of p at every degenerate p; 
in other words, EIG{p^^H,/3) is, like Pe,s, unstable against small perturbations of the 
Hamiltonian. (And, as with Pe,s, this fact, quite independently of the considerations 
on behalf of GAP-measures in Sections HI and suggests against using EIG{pp) as a 
thermal equilibrium distribution.) Moreover, EIG{p) is highly concentrated, generically 
on a one-dimensional subset of ^{Jif), and in the case of a finite-dimensional Hilbert 
space Jif fails to be absolutely continuous relative to the uniform distribution u,y{^jf) on 
the unit sphere. 

For further discussion of families p{p) of measures satisfying the analogues of Prop- 
erties 1-3, see Section ini21 

6.1.2 An Extremal Distribution 

Here is another distribution on associated with the density matrix p. Let the random 
vector ^ be 

VI/ = Vp^P' (42) 

pgspoc(p) 

the '^p being independent random vectors with distributions In case all eigen- 

values are nondegenerate, this means the coefficients of \E', \E' = Ylm'^nW), have 
independent uniform phases but fixed moduli \Zn\ = ^Jp^ — in sharp contrast with the 
moduli when is G'AP(p)-distributed. And in contrast to the measure EIG{p) consid- 
ered in the previous subsection, the weights pn in the density matrix now come from the 

^The relevant condition for (c) follows from nondegeneracy but is weaker: it is that the eigenvalues of 
pi and p2 are multiplicatively independent, in the sense that Pn^Pn] — Pm\pm\ can occur only trivially, 
i.e., when pn^ = pm\ and pl^^ — Pml- In particular, the nondegeneracy of pi and p2 is irrelevant. 

^ A property weaker than (c) does hold for EIG{p) also in the case of the degeneracy oi p — piiS)p2- if 
the orthonormal basis {|?2)} used in the definition of conditional wave function consists of eigenvectors 
of p2, then the distribution of the conditional wave function is EIG{pi). 
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fixed size of the coefficients of \1/ when it is decomposed into the eigenvectors of p, rather 
than from the probabihty with which these eigenvectors are chosen. This measure, too, 
is stationary under any unitary evolution that leaves p invariant. In particular, it is 
stationary in the thermal case p = p,X''^H,pi and for generic H it is an extremal station- 
ary measure as characterized in Section 13.2.11 in fact it is, in the notation of the last 
paragraph of Section 13.2.11 Ap with p{En) = (1/2') exp(— 

This measure, too, is highly concentrated: For a Hilbert space J€' of finite dimension 
k, it is supported by a submanifold of real dimension 2k — m where m is the number of 
distinct eigenvalues of H, hence generically it is supported by a submanifold of just half 
the dimension of J^. 

6.1.3 The Distribution of Guerra and Loffredo 

In PO], Guerra and Loffredo consider the canonical density matrix pis for the one- 
dimensional harmonic oscillator and want to associate with it a diffusion process on 
the real line, using stochastic mechanics [1311]. Since stochastic mechanics associates a 
process with every wave function, they achieve this by finding a measure pp on ^(L^(]R)) 
whose density matrix is p^. 

They propose the following measure Pj3, supported by coherent states. With every 
point {q, p) in the classical phase space of the harmonic oscillator there is associated 
a coherent state 

^,,p(x) = {2naY'/' exp (-^^^ + ^xp - ^pq^ (43) 

with 0"^ = h/2mLJ, thus defining a mapping C : ^ =5^(L^(]R)), C{q,p) = ipg^p. Let 
H{q,p) = p^ /2m + ^muj'^q^ be the classical Hamiltonian function, and consider the 
classical canonical distribution at inverse temperature /?', 

p^^ridq X dp) = ^e-'^'^^'^^PUqdp, Z' = [ dqdpe'^'^^'^^'P^ . (44) 

Let (3' = . Then pp = C^p'^?^^ is the distribution on coherent states arising from 

^ciass_ rjj^g density matrix of Pjs is pp [TU] . 

This measure is concentrated on a 2-dimensional submanifold of t5^(L^(M)), namely 
on the set of coherent states (the image of C) . Note also that not every density matrix p 
on L^(]R) can arise as the density matrix of a distribution on the set of coherent states; 
for example, a pure state p = can arise in this way if and only if -0 is a coherent 

state. 
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6.1.4 The Distribution Maximizing an Entropy Functional 

In a similar spirit, one may consider, on a finite-dimensional Hilbert space the 
distribution 'j{dtp) = f{jp)uy(^^){d'ip) that maximizes the Gibbs entropy functional 

^[f] = - j u{d^P) /(^)log/(^) (45) 
under the constraints that 7 be a probability distribution with mean and covariance 

PM',H,i3- 

/ > (46a) 
m(#) fiip) = 1 (46b) 

u{dij) f{tp) \^) = (46c) 



u{dtlj)f{tlj)\^){^\ = p,^,H,f3. (46d) 



A standard calculation using Lagrange multipliers leads to 

/(^) = exp(^|L|^) (47) 



with L a self-adjoint matrix determined by (j46b|l and (j46dj) : comparison with (jlSbjl 
shows that 7 is not a GAP measure. (We remark, however, that another Gibbs entropy 
functional, ^'[/] = — j^jX{dip) f^ip) \ogf{ip), based on the Lebesgue measure A on 

instead of u^'(jf), is maximized, under the constraints that the mean be and 
the covariance be p, by the Gaussian measure, f{jp)\{d%l)) = G{p){dip).) There is no 
apparent reason why the family of 7 measures should be hereditary. 

The situation is different for the microcanonical ensemble: here, the distribution 
ue,5 = GAP{pe,s) that we propose is in fact the maximizer of the appropriate Gibbs 
entropy functional Which functional is that? Since any measure jidip) on ^{J^) 
whose covariance matrix is the projection pE,5 = const. l[E,E+s\iH) must be concentrated 
on the subspace J^e,5 and thus cannot be absolutely continuous (possess a density) 
relative to we consider instead its density relative to = m_b,5, that is, we 

consider ^{difj) = f{jp)uE,5{d'ip) and set 



^"[/] = - j ns,5(#)/(^)log/(^). (4J 
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Under the constraints that the probabihty measure 7 have mean and covariance pE,5^ 
^"[/] is maximized by / = 1, or 7 = ue^s'-, in fact even without the constraints on 7, 
^"[/] is maximized by / = 1. 

6.1.5 The Distribution of Brody and Hughston 

Brody and Hughston [5] have proposed the following distribution n to describe thermal 
equilibrium. They observe that the projective space arising from a finite-dimensional 
Hilbert space, endowed with the dynamics arising from the unitary dynamics on Hilbert 
space, can be regarded as a classical Hamiltonian system with Hamiltonian function 
H{Cijj) = {iplHlip) / (tpltp) (and symplectic form arising from the Hilbert space struc- 
ture). They then define n to be the classical canonical distribution of this Hamiltonian 
system, i.e., to have density proportional to exp(— /?if(C'?/')) relative to the uniform vol- 
ume measure on the projective space (which can be obtained from the symplectic form 
or, alternatively, from M.j^(^) by projection from the sphere to the projective space). 
However, this distribution leads to a density matrix, different from the usual one pp 
given by (jHI), that does not describe the canonical ensemble. 

6.2 A Uniqueness Result for GAP{p) 

As EIG{p) is a family of measures satisfying Properties 1-3 for most density matrices p, 
the question arises whether there is any family of measures, besides GAP{p), satisfying 
these properties for all density matrices. We expect that the answer is no, and formulate 
the following uniqueness conjecture: Given, for every Hilbert space Jif and every density 
matrix p on Jif, a probability measure p{p) on S^{J^) such that Properties 1-3 remain 
true when GAP{p) is replaced by p{p), then p{p) = GAP{p). In other words, we 
conjecture that p = GAP{p) is the only hereditary covariant inverse of (0). 

This is in fact true when we assume in addition that the mapping p : p p{p) is 
suitably continuous. Here is the argument: When p is a multiple of a projection, p = 
(dim^')~^P^^/ for a subspace Jif' C then p{p) must be, by covariance U^p{p) = 
p{UpU*), the uniform distribution on J>^{Jf"), and thus p{p) = GAP{p) in this case. 
Consider now a composite of a system (system 1) and a large heat bath (system 2) with 
Hilbert space = M'\ ® and Hamiltonian H = Hi ® I2 + Ii ® and consider 
the microcanonical density matrix pe^s for this system. By equivalence of ensembles, we 
have for suitable /5 > that pe,s ~ P.%',H,f3 = Pp'^ ® pf^ where p^*^ = p.%i,Hi,f3- By the 
continuity of p and GAP, 

Kp? ® P?) ^ PipE,s) = GAP{pE,5) ^ GAP{p^^^ ® pf) . 
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Now consider, for a wave function \1/ with distribution ^-^Pj^ ®pf ) respectively GAP{p^ ® 
pf^), the distribution of the conditional wave function \E'i: by heredity, this is p{p^p^) 
respectively GAP{p^^^). Since the distribution of \E'i is a continuous function of the 
distribution of \E', we thus have that n{p^p^) ~ GAP{p^p^). Since we can make the degree 
of approximation arbitrarily good by making the heat bath sufficiently large, we must 
have that p{p^p^) = GAP{p^^^). For any density matrix p on J^i that does not have zero 
among its eigenvalues, there is an Hi such that p = p^^^ = exp{—(3Hi) for (3=1, 
and thus we have that p{p) = GAP{p) for such a p; since these are dense, we have that 
/i(p) = GAP{p) for all density matrices p on Since J^i is arbitrary we are done. 

6.3 Dynamics of the Conditional Wave Function 

Markov processes in Hilbert space have long been considered (see |lj for an overview), 
particularly diffusion processes and piecewise deterministic (jump) processes. This is 
often done for the purpose of numerical simulation of a master equation for the density 
matrix, or as a model of continuous measurement or of spontaneous wave function 
collapse. Such processes could arise as follows. 

Since the conditional wave function arises from the wave function (52 1 V^) by in- 
serting a random coordinate Q2 for the second variable (and normalizing), any dynamics 
(i.e., time evolution) for Q2, described by a curve t ^— Q2(t) and preserving the quantum 
probability distribution of Q2, for example, as given by Bohmian mechanics jB], gives 
rise to a dynamics for the conditional wave function, t 1— > \E'i(t) = A/'(t) {Q2{t)\ip(t)) , 
where ip{t) evolves according to Schrodinger's equation and A/'(t) = || (<52(^)|V^(^)) ||^^ is 
the normalizing factor. In this way one obtains a stochastic process (a random path) 
in y{J^i). In the case considered in Section IH in which ^ corresponds to a large 
heat bath, this process must have GAP{p,^{^Hi,i3) as an invariant measure. It would be 
interesting to know whether this process is approximately a simple process in S^{Jifi), 
perhaps a diffusion process, perhaps one of the Markov processes on Hilbert space con- 
sidered already in the literature. 

7 The Two-Level System as a Simple Example 

In this last section, we consider a two-level system, with = and 

H = Ei\l){l\ + E2\2){2l (49) 

and calculate the joint distribution of the energy coefficients Zi = and Z2 = (2|\I') 

for a G'AP(p/3)-distributed \E' as explicitly as possible. We begin with a general finite- 
dimensional system, = C^, and specialize to k = 2 later. 
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One way of describing the distribution of \l/ is to give its density relative to the 
hypersurface area measure u on y{C^); this we did in (HHD- Another way of describing 
the joint distribution of the Z„ is to describe the joint distribution of their moduh |Z„|, 
or of l^nP; as the phases of the Z„ are independent (of each other and of the moduh) 
and uniformly distributed, see (j221)- 

Before we determine the distribution of |^„p, we repeat that its expectation can be 
computed easily. In fact, for any (p G Jif we have 



Thus, for = \n),we obtain E|Z„,p = e"'^^"/tr e"'^^. 

For greater clarity, from now on we write Z^^^ instead of A relation similar 
to that between GAP{p), GA{p), and G{p) holds between the joint distributions of the 
IZ^^^I"^, of the iZ^'^p, and of the IZ^]"^. The joint distribution of the IZ^]"^ is very sim- 
ple: they are independent and exponentially distributed with means p„ = e~^^"/Z(/?). 
Since the density of GA relative to G, dGA/dG = Yin l^nl"^, is a function of the moduli 
alone, and since, according to (j221), GA = GAphases x GA^oduH, "we have that 

C'^^moduli ^ ^ I \ G^moduli • 
n 

Thus, 

P(|Zf^P G ds„ e ds,) = '^ + --- + "^ exp(- \2-)ds,--- ds,, (50) 

where each Sn G (0, oo). Finally, the \Z^^^\'^ arise by normalization, 

I 7GA\2 

\vGAP\2 _ \^n I /p-.N 
n' " 

We now specialize to the two-level system, k = 2. Since IZ^^^I"^ + \Z2^^\^ = 1, it 
suffices to determine the distribution of \Z^^^\'^, for which we give an explicit formula 
in ()52c|l below. We want to obtain the marginal distribution of (jKT|) from the joint 
distribution of the l^^'^P in (0, oo)^, the first quadrant of the plane, as given by (j^UI) . 
To this end, we introduce new coordinates in the first quadrant: 

_ ^1 \ _ , 

S — ; , A — Si -\- S2, 

Si + S2 

where A > and < s < 1. Conversely, we have Si = sX and S2 = (1 — s)A, and the 
area element transforms according to 

d{s,X) 



dsi ds2 



dsdX = X ds dX. 
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Figure 1: Plot of the distribution density function /(s) of defined in fl52cp . for 

various values of the parameter 6 = exp(/3(£'2 — -Ei)): (a) 6 = 1/3, (b) S = 1/2, (c) 6=1, 
(d) 6 = 2, (e) 6 = 3. 

Therefore, using 

oo 


we obtain 

oo 

p(|Zf^^|2 e rfs) = rfs j rfA^A2exp(-A(e^^is + e'3^^(l-s))) = (52a) 



= ^-^{e'''s + e^^^{l-s))-'ds= (52b) 

= [ais + a2{l - s))~^ ds =: f (s) ds , < s < 1, (52c) 

with tti = {6-\6-^ + l)/2)i/3 and aa = (<5(5 + l)/2)i/3 for 6 = exp{(3{E2 - ^i)). The 
density / of the distribution (j52cp of \Z'^'^^\'^ is depicted in Figure 1 for various values 
of 6. For 5 = 1, / is identically 1. For 5 > 1, we have a2 = 6ai > ai, so that 
ais + 02(1 — s) is decreasing monotonically from 0:2 at s = to ai at s = 1; hence, / is 
increasing monotonically from to For 5 < 1, we have 0:2 < ai, and hence / is 
decreasing monotonically from to aj"'^. In all cases / is convex since /" > 0. 
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